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ABSTRACT 

Optical CCD imaging with Ha and [SII] filters and spectroscopic observa- 
tions of the galactic supernova remnant G85.9-0.6 have been performed for the 
first time. The CCD image data are taken with the 1.5m Russian-Turkish Tele- 
scope (RTT150) at TUBiTAK National Observatory (TUG) and spectral data 
are taken with the Bok 2.3 m telescope on Kitt Peak, AZ. 

The images are taken with narrow-band interference filters Ha, [SH] and 
their continuum. [SH]/Ha ratio image is performed. The ratio obtained from 
[SII]/Ha is found to be ~0.42, indicating that the remnant interacts with HII 
regions. G85. 9-0. 6 shows diffuse-shell morphology. [SH]AA6716/6731 average fiux 
ratio is calculated from the spectra, and the electron density is obtained to be 
395 cmT^. Prom [0IIl]/H/3 ratio, shock velocity has been estimated, pre-shock 
density of ric = 14 cm~^, explosion energy of £^ = 9.2 x 10^° ergs, interstellar 
extinction of E{B — V) — 0.28, and neutral hydrogen column density of N{HI) — 
1.53 X 10^^ cm~^ are reported. 

Subject headings: G85.9-0.6; optical observation; CCD image; spectra 
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Introduction 



Galactic supernova remnants (SNRs) are now as many as 274 flGreen II2009I ). Most 
of tliem are observed in radio wavelengths from their non-thermal synchrotron emission 
despite their very few optical observations. The reason for this is because of the dust 
and gas around the galactic disk stopping the optical light coming from them. However, 
only recently optical observations of SNRs can be performed by using narrow-ban c 



interference filters and providing rather long posing intervals (e.g. 



Mavromatakis et al. 



2001 



Mavromatakis 



2003 



Boumis et al. 



Fesen et a. 



2002, 



2009 



1997, 



2008; 



GoketaL 



20081 ). Recently, X- ray satellites li 



in X-ray band (e.g. 



Bamba et al. 



^e Chandra, XMM-Newton , Suzaku observed some SNRs 



200ol : 



Katsuda et al. 



20091). 



The galactic supernova remnant G85.9-0.6 (a = 20^58'"40^ ^ = 4 4°53') was discovered 



in the Canadian Galactic Plane Survey (CGPS) (IKothes et al. 
radio surface brightness at 1 GHz of 2 x 10^^^ Wattm~'^Hz^^sr~^. 



2001). G85.9-0.6 has a 



Kothes et al. 



(1200lh 



has detected no HI features corresponding to the SNR. This and the low radio and X-ray 
brightness, suggest expansion in a low-density medium. The SNR may lie in the low-density 
region between the local and Perseus spiral arms, at a distance of about 5 kpc. The most 
likely event that would produce an SN R in such a low dens ity medium might be a Type 
la supernova. It has a radius of ~ 12' ( Kothes et al. II2OOII ). G85.9-0.6 has a diffuse type 



structure in X-ra y region, where its e 



order of 1.0 Mq (IJackson et al. 



20081). 



ectron density is 0.20 cm ^ and the ejecta mass is at 



[SIIj/Hg ratio is an iraportant criter ia to distinguish HII regions and SNRs 



( iMathewson &: Clarke 



1972 



Dopita 



19971 ). In HII regions most sulphur atoms are in the 



form of SHI, since there is a strong photo-ion ization and [SII]/Ha ratio is expected to be 



in the r ange of ~ 0.1 — 0.3 (lOsterbrock 



1[ 

0.1-0.5 fiDopita et al. iboool : iKewlev et al. 



1989), h owever recent models give this range to be 



2001 



), for SNRs on the other hand, it is 0.5-1.0 



-A- 



flRavmond Ill979l : IShuU fc McKee Ill979f ) 



In this work, we selected the SNRs (G16.2-2.7, G17.8-2.6, G36.6+2.6, G83.0-0.3 and 
G85.9-0.6) that have no optical observation in literature so far and have proper coordinates 
to observe at TUG. We take their Ha images with 300 seconds exposure time. 

Theoretical models predict that, in SNRs, the temperature of the regions that give 
optical emission should be at ~ 10^ K, therefore, the intensity of Ha emission line should 
be larger compared to the other emission lines. Fo r this reason. Ha interference filters are 



preferred for determining SNRs (jOsterbrock 



19891 ). With Ha filter we found radiation only 



from G85.9-0.6 and studied its optical image and spectrum. The observations and results 
are presented in sections 2. and 3., respectively. In the last part we discuss the results. 



2. Observations 
2.1. Imaging 

The CCD images of G85.9-0.6 are obtained with the 1.5m RTT150 at TUG, Antalya, 
Tiirkiye. The data are taken with a 2048x2048 pixel CCD447, with a pixel size of 15/im 
resulting in a 13.5 arcmin x 13.5 arcmin field of view, with a plate scale of 0".39 pixel~^. 
Images are obtained with Ha, [SII] and their continuum filters within 900 seconds exposure 
times. Continuum filters are used to eliminate contamination from background starlight. 
Bias frames and dome fiats are also observed for data reduction. The observations and 
interference filter characteristics are given in Table 1. The coordinates quoted in this work 
refer to epoch 2000. The image reduction is performed by using standard IRAF (Image 
Reduction Analysis Facility) routines. 

Backgrounded, fiat-fielded, trimmed and continuum subtracted [SII] and Ha images 
are prepared. In Figure 1, Ha, Ha cont., Ha-Ha cont. and [SII], [SII] cont., [SII] -[SII] cont. 
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images are presented. The images have been smoothed to suppress the remaining of the 
imperfect continuum subtraction. 



2.2. Spectroscopy 

An optical spectrum of SNR G85.9-0.6, covering 3750 - 7250 A range, was obtained on 
UT 2009 June 16 at the Bok 2.3 m telescope on Kitt Peak, AZ, using the Boiler and Chivens 
Spectrograph at the Richey- Chretien f/ 9 focus together with a thinned back-illuminated 
1200x800, 15/im pixel Loral CCD. We used a 1st order 400 g/mm grating blazed at 4889 A. 
With a slit width of 2.5 arcseconds a typical resolution of 6 A was achieved at a reciprocal 
dispersion of 2.78 A / pix on the CCD. Spectrums are taken on the image from three 
different intensive emission regions far from the stars. The coordinates of these regions and 
exposure times are listed in Table 2. 

The data reduction is carried out by using the standard IRAF routines. Each frame 
is bias subtracted and for flat- field correction halogen lamps are used. Then with apall 
package under noao.twodspec.apextract, each spectrum is cleaned from external factors like 
sky and cosmic rays. The apall package selects an aperture for the target spectrum, and 
two apertures for the sky spectrum. The apertures for the sky spectrum are selected a 
dozen or so pixels, not illuminated by the object, above and below the target spectrum. 
Apall takes the average of the two sky spectrum and subtracts it from the target spectrum 
giving a sky subtracted target spectrum. Identify, reidentify, hedit and dispcor packages 
under noao . onedspec are used for wavelength calibration. For flux calibration, sensfunc 
package under noao. oned s pec an d calibrate package under noao. twodspec. long slit are used 



as described by 



Massey 



( 1l997l ). The spectra is dereddened by using dereddened task 



under noao. onedspec. The spectrophotonaetric standard star HR5501 is observed for flux 



calibration of the spectra (IHamuy et al 



19921 ). For wavelength calibration, however, 
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helium-neon (He-Ne) arc lamps are used. 

The spectra obtained for three region are shown in Figure 2. In the optical spectra 
of SNRs; Ha, Rp and [OIII]A4959, A5007, [NII]A6548, A6584, [SII]A6716, A6731 known as 
forbidden lines, are expected. Taking the temperature 10^ K which is convenient for optical 
radiation, the density of the medium, the energy of the supernova explosion, interstellar 
absorption, neutral Hydrogen column density, shock wave velocities and all other related 
physical parameters are derived from these lines. 



3. Results 

In Figure 3, we present [SII]-[SII]cont./Ha-Ha; cont. image of the remnant G85.9-0.6. 
Continuum substraction and image division are done with IRAF imarith package. Imarith 
package subtracts and divides the number of photons in each pixel of the image obtained 
with the above mentioned filters. Calculated [SII]/Ha ratios and their errors in parenthesis 
taken from 10 different regions on the image are listed in Table 3, the average of these ratios 
is 0.42±0.04. 

Figure 4 presents the correlation between Ha emission and the radio emission of G85.9- 



0.6 at 1.4 GHz (radio contours are obtained from http: //skyview.gsfc. nasa.gov/cgi-bin/query. pi ) 



The western diffuse emission is well correlated with the radio emission. 

The flux values of the spectra shown in Figure 2 are obtained by normalizing them with 
respect to Ha flux value. With these flux values and taking the electron temperature at 
~ 10^ K, some basic parameters of SNR G85.9-0.6 are calculated and listed in Table 4. The 
fluxes corrected for interstellar extinction are normalized to I(Ha) = 100. For three region, 
[SII]/Ha, [SII]AA6716/6731, [0III]/H/3 ratios and the parameters that are obtained from 
these ratios found to be comparable. The average values are considered while discussing the 
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results. 



To find t hese paramet 



is calculated (lOsterbrock 



ers fi rst, through [S11]AA6716/6731 flux ratio the electron density 



19891). Shock velocity is estimatec 



described by (I Cox fc Raymond 



1985 



Raymond et al. 



1988 



Tom the ratio o 


F [Oil] 


Hartipan et al. 


1987) 



I]/H/? as 



Assuming that the remnant is still in the adiabatic phase of its evolution, the pre-shock 
cloud density jir can be calculated by u sing the relationship n[SII] = 45nc {Vs/IOOY cm~^ 



( iDopita 



1979 : 



Fesen &: Kirshner 



19801 ). Where n[Sll] is the density calculated using the 
[SII] emission line from the optical spectra. This equation is valid if the transverse magnetic 
field is sufficiently low suc h that the magnetic pr e ssure does not dominate the post-shock 



pressure as is discussed by iDopita fc Sutherland 



fll996h . 



e.g. 



Shock wave energy E i s given by the equation 2 x 10^^ [3 ^ ric (V^s/100)^r^ ergs (see 



McKee fc Cowie 



19751 ) where is the radius of the shock wave. The factor (3 is 



approximately equal to 1 at the blast wave shock. 

The interstellar extinction {E{B — V")), absorption {Av), and neutral hydrogen 
column density {N{HI)) are calculated through the relati ons E(B — V ) = 0.664c, 



Av = 3.1 xE(B-V) and 



(jPredehl fc Schmitt 



N{HI) = 5.4 X 10^^ X E{B - V) flKaler 



19761 ). flAUer 



19841 ) and 



19951 1. Where c is 1/0.331 log[{Ha/ H (3) /?,]. 



Furthermore, [SII] /Ha ratio is also calculated from the spectrum. In this calculation, 
A6716 and A6731 doublet line fiuxes of [SII] is used. This ratio, for three coordinates whose 
a and (3 are given in the Table 2, is calculated at its average value to 0.46. 
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Discussion and Conclusions 



In this work, the imaging and spectroscopic observations of the central part (13' x 13 ) 
of the galactic supernova remnant G85.9-0.6 is observed for the first time in optical band. 
The images show diffuse emission structure. As seen from the Figures 1,3, the dominant 
emission extends from south-west to north. The diffuse emission centered at a = 20'^59™00^ 
5 = 44°52'50" that covers an area 6 x 8 is very well correlated with radio emission (see 
Fig.4). 

As seen from the Tablo 3, [SII]/Ha ratios take values between 0.2-0.6. This ratio is 
~0.5, 0.6 on the places where dominant diffuse emission comes indicating ionization of 
the shock heated gas resulted from collision. On the other hand, the ratio of ~0.2-0.3 
on the places where emission is faint indicating that resulted emission originate from 
photoionization mechanism. [SII]/H« ratio obtained from imaging is at 0.42 which is in 
agrement with spectral measurements of 0.46. These values imply that diffuse emission of 
G85.9-0.6 may be associated with HII region. 

From the spectra, we see the presence of Ha, H/3, [OIII]A4959, A5007, [NII]A6548, A6584 
and [SII]A6716, A6731, in addition a weaker [OI]A6300 and [OIII]A4363 emission lines. 
The presence of [OI]A6300 line emission indicates that the emission is originating 
from shock heated gas. Theoretical models predict electron temperature through 
[OIII]J(A4959 + 5007/A4363)ratio. Since [OIII]A4363 emission hue is v ery close to Hg I 



19891 ). Therefore, the 



A4358, it is almost impossible to differentiate them (jOsterbrock 
electron temperature (T) assumed to be 10^ K in this work. 

We detected a weak [OIII]/ H/? emission line ratio at 1.18 i ndicating an oxygen- de ficient 



remnant. Theoretical models of 



Cox fc Raymond 



fll985h and 



Hartigan et al. 



fll987h give 



that for shocks with complete recombination zones this value is ' ;^6, while this limit is 



exceeded in case of shocks with incomplete recombination zones (iRaymond et al. 



19881 ). 
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Our values measured from spectra indicates that complete recombination takes place 
at those regions. Therefore, considering our measurements and the theoretical models 
mentioned above we estimated the velocity of the shock to be 80 kms~^. 

Since the G85.9-0.6 lie s in the region betw een the local and Perseus spiral arms, it has a 



reliable distance, of 5 kpc (IKothes et al. 



20011 ) ■ For this distance and 9 = 24' we calculated 



shock radius (r^) to be ~ 17.4 pc. Explosion energy and hydrogen column density of the 
G85.9-0.9 for this radius and estimated shock velocity is calculated, to be ~ 9.2 x 10^° ergs 
and ~ 1.53 x 10^^ cm^^ respectively. This much energy can be considered as the typical 
energy released from a supernova explosion. 

Typical value of t 



1.5 ( Frank et al. 



;he sulfur line ratio [SII](A6716/A6731) for SNRs should be close to 
2002). This ratio is calculated to be 1.13 indicating low electron density 
which is calculated as 395 cm~^ for this region. 

We note here that, providing rather long exposure times and examining its surrounding 
region will give more information about the morphology and angular size of the G85.9-0.6 
in optical band. The images taken with different interference filters like [Nil], [OIII] will 
give more information about the inhomogeneities and density variations of the medium. 
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Fig. 1. — Ha, Ha cont., Ha-Ha cont., [SH], [SH] cont. 
G85.9-0.6. Each field of view is 13.5 arcmin x 13.5 arcmin 



and [SH]-[SH] cont. images of 
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Fig. 2. — The long slit spectra of G85.9-0.6 for area 1 (a), area 2 (b) and area 3 (c) 
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Fig. 3. — The negative of the [SII]/Hq; image of G85.9-0.6. The image has been smoothed 
to suppress the residuals from the imperfect continuum subtraction. Field of view is 13.5 
arcmin x 13.5 arcmin 
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Fig. 4. — NVSS 1.4-GHz VLA radio continuum map superposed on the continuum subtracted 
and smoothed Ha image of G85.9-0.6 
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Table 1: The journal of observations. 

Filter Wavelength Obs. Date (UT) Exp. times (s) 

(FWHM) (A) 

Ha 6563 (88) 01 June 2009 900 

Ra cont. 6446 (130) 01 June 2009 900 

SIX 6728 (70) 01 June 2009 900 

SIX cont. 6964 (300) 01 June 2009 900 
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Table 2: The spectroscopic observations. 



Area Slit Center (J2000) Exposure times (s) 

1 a = 20'*59'"04^ 600 
5 = 44^52' 19" 

2 a = 20^58^49^ 600 
S = 44055'32" 

3 a = 20'*58"^32^ 600 
5 = 44058'24" 
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Table 3: [SII]/Ha ratios and their errors for different regions of tlie G85.9-0.6. 



A 

Area 


coordinate J (2000) 


ri~1TTl /TT 

[SIIJ/Hq! 




n 


. () 




1 


20'*58"*29" 


, 44°58'37" 


0.6 (± 0.1) 


2 


20/^58^335 


, 44058'30" 


0.5 (± 0.1) 


3 


20''58"'49" 


, 44055'07" 


0.5 (± 0.1) 


4 


20'*58"'53" 


, 44°53'27" 


0.5 (± 0.1) 


5 


20^58^^48" 


, 44055'20" 


0.4 (± 0.1) 


6 


20'*58'"29' 


, 44°48'20" 


0.3 (± 0.1) 


7 


20''58"'18" 


, 44°56'02" 


0.3 (± 0.1) 


8 


20''58"'11" 


, 4405l'l7" 


0.2 (± 0.1) 


9 


20^59™03" 


, 4405/29" 


0.3 (± 0.1) 


10 


20''59"'04^ 


, 44°53'53" 


0.6 (± 0.1) 



Table 4: This table gives the relative line fluxes and the parameters obtained from spectra for three different areas of the 
G85.9-0.6. (F) shows fluxes uncorrected for interstellar extinction, (I) shows fluxes corrected for interstellar extinction. 
Numbers in parentheses represent the signal-to-noise ratio of the quoted values. All fluxes normahzed F{Ha)— 100 
and l{Ha)~100. The errors of the emission line ratios and other parameters are calculated through standard error 
calculation. 
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100 T^HrvV 


-1 nn^ 
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V cLl Uco 








ArGci 1 




Area 2 




A rt^a -< 
-Til cd O 














"p 

r 




F 


I(S/N) 


r 










A rt^c 9 
CcL Z 


rvl Cd-O 


4861 


28 


52(58) 


25 


48(14) 


20 




34(26) 


[SII] /Ha 


0.44±0.01 


0.45±0.01 


0.50±0.09 


4959 [OIII] 


8 


15(40) 


11 


21(7) 


8 




13(9) 


[SII]A6716/6731 


1.09±0.08 


1.14±0.03 


1.17±0.2 


5007 [OIII] 


24 


44(54) 


19 


35(10) 


18 




29(23) 


[0III]/H/3 


1.14±0.14 


1.16±0.20 


1.23±0.16 


6300 [01] 


7 


17(28) 


9 


13(4) 


5 




6(4) 


Ne (cm~^) 


470±120 


380±40 


336±26 


6548 [Nil] 


10 


10(28) 


10 


11(5) 


13 




14(11) 


lie (cm~^) 


16±3 


13±1 


12±1 


6563 Ha 


100 


100(245) 


100 


100 (30) 


100 




100(124) 


c 


0.23±0.02 


0.38±0.04 


0.67±0.07 


6584[NII] 


29 


33(56) 


33 


31(9) 


33 




32(31) 


E(B-V) 


0.15±0.02 


0.25±0.03 


0.45±0.05 


6716 [SII] 


25 


23(40) 


24 


24(8) 


23 




27(7) 


Av 


0.47±0.05 


1.18±0.12 


1.40±0.14 


6731 [SII] 


20 


21(35) 


20 


21(7) 


19 




23(4) 











